Optimizing the regimes of Advanced LIGO gravitational wave detector for multiple 

source types 
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We develop here algorithms which allow to find regimes ol signal-recycled Fabry-Perot-Michelson 
interferometer (for example, Advanced LIGO), optimized concurrently for two (binary inspirals + 
bursts) and three (binary inspirals + bursts + millisecond pulsars) types of gravitational waves 
sources. We show that there exists a relatevely large area in the interferometer parameters space 
where the detector sensitivity to the first two kinds of sources differs only by a few percent from 
the maximal ones for each kind of source. In particular, there exists a specific regime where this 
difference is w 0.5% for both of them. Furthermore we show that even more multipurpose regimes 
are also possible, that provide significant sensitivity gain for millisecond pulsars with only minor 
sensitivity degradation for binary inspirals and bursts. 



I. INTRODUCTION 

Within the last decade we have witnessed a very signif- 
icant progress in experimental gravitational wave (GW) 
astronomy. All the ground-based interferometric GW an- 
tennae such as LIGO [1] in the USA, VIRGO [2] in Italy 
GEO600 [3] in Germany and TAMA300 [4] in Japan have 
been commissioned to operation and started to record sci- 
entific data. Nevertheless, no signs of gravitational waves 
were found thus far in this data which is, as we under- 
stand now, quite reasonable as it agrees with moderately 
optimistic predictions of the astrophysicists on the rate of 
measurable events within the limits of antennae detection 
range. This possibility was realized by GW community, 
and work on design of the next, second generation of GW 
antennae went on in parallel with efforts in enhancement 
of the first generation ones. A pioneer amongst the sec- 
ond generation GW detectors will become an American 
Advanced LIGO project whose construction should start 
in 2010 [5]. It is planned to have sensitivity more than an 
order of magnitude higher than its predecessor. Such a 
dramatic increase will be provided by significantly lower 
seismic noise level due to new active antiscismic isolation, 
use of higher quality optics and lower level of quantum 
noise. 

The main difference between the Initial LIGO and Ad- 
vanced LIGO designs that is crucial for lowering this 
noise is the use of signal recycling (SR) technique. Its im- 
plementation in contemporary detector setup is relatively 
easy as it requires to install only one additional mirror in 
the interferometer output port. This mirror reflects side- 
band signal field coming out of the interferometer back to 
the arm cavities or "recycles" it. However, the dynamics 
and quantum noise properties of the interferometer be- 
come much richer and thus provide more freedom in ad- 
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justment its sensitivity curve to fit the current research 
goals. In particular, it was stressed by A. Buonnano and 
Y. Chen [6] that the optical system composed of the SR 
cavity and the arm cavities forms a composite resonant 
cavity whose eigenfrequencies and quality factors can be 
controlled by the position and reflectivity of the SR mir- 
ror, thus increasing or decreasing the storage time of the 
signal inside the cavity. Moreover, as in Advanced LIGO 
it is planned to increase the optical power, circulating 
in the arm cavities approximately 80 times with respect 
to initial LIGO, an optomechanical interaction between 
laser field and mirrors will significantly influence the dy- 
namics of the test masses turning them from free bodies 
(within the detection frequency band: ~ 10 -f- 10 4 Hz) 
into oscillators with eigenfrequency falling into the de- 
tection band. This effect, known as ponderomotive rigid- 
ity [6, 7, 8, 9, 10] arises when off-resonant optical field 
creates an effective restoring force originating from radi- 
ation pressure, which occurs to be a function of mirror 
displacement. In this situation one can say that optical 
field creates a frequency dependent mechanical rigidity 
[6, 9, 11, 12, 13, 14]. 

As a result, the quantum noise spectral density of sig- 
nal recycled interferometers (SRI) can be tuned to pro- 
vide the best sensitivity for different gravitational wave 
sources. So far, it was supposed that in order to reach 
good sensitivity for each of the source types, totally dif- 
ferent strategies should be used, which correspond to dif- 
ferent optical parameters sets. Detection of gravitational 
waves from inspiraling neutron star binaries (NSNS) re- 
quires, for example, that the noise spectral density have 
to be as small as possible at low and medium frequencies, 
/ < 100 Hz. If one is interested in narrow-band detec- 
tion of GWs from the source with well defined center fre- 
quency /puis, such as high frequency pulsars, the optimal 
regime of the interferometer will be absolutely different: 
evidently, one should choose such set of optical param- 
eters that provides minimum to noise at / pu i s - Searches 
of GWs from supernovae bursts or stochastic relic grav- 
itational radiation which analytical waveforms are not 
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FIG. 1: Relative improvement in Signal-to-Noise-Ratio 
Pburst/Pburst f° r bursts of GW radiation as a function of rel- 
ative deterioration in SNR for neutron stars (NS) binaries 
1 - Pns/Pns x - Here p®*% is the value of SNR for GW burst 
sources when interferometer is tuned for reaching maximum 
SNR for NS binaries. 
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FIG. 2: Relative improvement in Signal-to-Noise-Ratio 
Ppuis/Pp^if for one of the high frequency pulsars (J0034-0534 
with rotational frequency fo ~ 532, 7 Hz) as a function of 
relative deterioration in SNR for neutron stars (NS) bina- 
ries 1 - pns/pns x - Here p®£f is the value of SNR for GWs 
from specified pulsar when interferometer is tuned for reach- 
ing maximum SNR for NS binaries. 



known require flat broadband noise curves [15, 16, 17]. 

Extensive studies of optimal modes of operation of 
Advanced LIGO interferometer most suitable for differ- 
ent specific kinds of sources of GWs were carried out 
in AdvLIGO Lab in extensive details [18]. However, 
in this paper we focus on finding multipurpose regimes 



of signal recycled interferometers which might provide 
good (sub-optimal) sensitivity simultaneously for differ- 
ent gravitational wave source types. The possibility to do 
this originates from the fact that classical noise budget 
for currently operating and future GW detectors masks 
low and medium frequency features of the quantum noise 
(see Fig. 5) that are mostly susceptible to variation of 
optical parameters of the interferometer. This fact, al- 
though being quite unpleasant for tunability of GW de- 
tector for specific sources is surprisingly advantageous 
for tuning the antenna to have high enough sensitivity to 
GWs from various types of astrophysical sources simul- 
taneously. As will be shown in subsequent sections the 
sensitivity changes relatively slowly within quite a wide 
range of main SRI optical parameters for different types 
of signals (GWs from inspiraling compact binaries, GW 
bursts, high-frequency pulsars etc.). And these areas for 
different sources significantly overlap that allows to find 
quasi-optimal regime for two or even three different GW 
sources simultaneously. It is shown by the example of Ad- 
vanced LIGO SRI that rather significant improvement in 
sensitivity to GW bursts and GWs from high frequency 
sources is possible at the cost of quite moderate dete- 
rioration of signal strength for compact binary systems. 
This is illustrated in Figs. 1 and 2 where the relative im- 
provement in Signal-to-Noise-Ratio (SNR) for GW bursts 
and high frequency pulsars, correspondingly, are plotted 
with respect to relative deterioration of SNR for NSNS, 
provided that one diverts optical parameters of SRI from 
the optimal ones for NSNS. It should be also emphasized 
here that in spite Advanced LIGO is used as an example 
for which we perform calculations, the results we obtain 
are general and applicable to all SRI limited by classical 
noise at low and medium frequencies. It is also instruc- 
tive to mention that our optimization includes only the 
most basic parameters and the results can be considered 
only as some preliminary guidelines for designing future 
generation of GW interferometers, while for more specific 
optimization of real device much larger parameter space 
should be considered and optimization over parameters 
that should be set before the device is built and the ones, 
that can be easily tuned in already operating detector 
should be performed separately. The above issues will be 
considered in future works. 

The paper is organized as follows. In Section II the 
brief consideration of quantum noise of signal-recycled 
Fabry-Pcrot-Michclson interferometer is performed. In 
Section III expressions for Signal-to-Noise Ratio and de- 
tection range for the gravitational-wave radiation from 
the inspiraling binary system are given and the numer- 
ical optimization procedure with respect to interferom- 
eter optical parameters is described. The quantitative 
and qualitative analysis of the obtained results of opti- 
mization against GWs from neutron star binaries is per- 
formed. In Section IV the sensitivity of interferometer to 
GW bursts is analyzed and regimes for simultaneous de- 
tection of these two types of sources are investigated. In 
section V this analysis extended also to high-frequency 
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Description 



Test bodies reduced mass 
Speed of light 
SRI arms length 
Laser frequency 
Mechanical frequency 
Circulating optical power 
Effective SRI half-bandwidth 
Effective SRI detuning 

Renormalized optical power 

Homodyne angle 

Total readout quantum efficiency (incl. losses) 



TABLE I: Notations used for characterizing quantum noise of 
SRI 
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FIG. 3: Principle optical scheme of signal-recycled interfer- 
ometer of Advanced LIGO GW detector 



quasi-periodic sources (pulsars). In Section VI the ob- 
tained results are discussed and some concluding remarks 
are given. Some notations and values of parameters used 
through this paper are listed in Table I. 



II. QUANTUM NOISE OF SIGNAL RECYCLED 
INTERFEROMETERS 



In Fig. 3 the schematic drawing of a signal recycled in- 
terferometer is presented. Here the additional signal re- 
cycling mirror (SRM) forms, together with the input test 
masses (ITMs) of arm cavities, an additional SR cavity 
which properties are defined by two parameters of SRM, 



namely its amplitude reflectivity p and detuning phase 
0src = [upl/c] mod 2ir gained by carrier light travelling 
one-way in the SR cavity (I is the length of SR cavity) . 

As demonstrated by A. Buonnano and Y. Chen [12], 
there exists one-to-one transformation ( "scaling law" ) be- 
tween the parameters of a single detuned Fabry-Perot 
(FP) cavity with one movable mirror and ones of SRI 
that allows to describe the optical behavior of it in terms 
of much more simple equivalent system such as FP cav- 
ity. According to "scaling law" for any SRI there exists 
a unique FP cavity with bandwidth 7 and detuning 6 
defined by formulae: 



(1 - p 2 )7AR.M 

1 + 2/3 cos 2</>src + p 2 

2/Q7ARM sin 2(/> SRC 
1 + 2/3 cos 20SRC + P 2 



(la) 
(lb) 



where 7arm = cT/AL is the half-bandwidth of arm FP 
cavities, that has the same optomechanical features and 
therefore the same sensitivity as the initial SRI. The effec- 
tive optical power circulating in the equivalent FP cavity 
should be twice as large as real optical power W circulat- 
ing in a single arm cavity. The same is referred to masses 
of input (ITM) and end (ETM) test masses of effective 
cavity: M eS = 2M. 

Below we will use also extensively the following conve- 
nient parameters: generalized bandwidth 



r / 2 1 x2 1 ~ 2pcos20 S RC + P 2 ,„ » 

r = V 7 + 2 = jARM\ h——: rr ; , ( 2a ) 

y 1 + 2pcos2^ S Rc + P 2 



and detuning phase 
8 



(3 = arctan ■ 



arctan 



1 



■2p 



1-p- 



sin 20 SRC 



(2b) 



The above expressions can be also easily reverted to ob- 
tain the SR cavity parameters: 



' 7ARM ~ 2 7ARM-T COS + r 2 

7arm + 27arm^ cos 13 + r 2 ' 



(3a) 



»SRC 



1-P 

arcsm [ tan p 

2p 

fl-p 2 
tt — arcsin [ — tan (3 



r < 7ARM , 
r > 7ARM • 

(3b) 



Using the "scaling law" approach, consider FP cav- 
ity with movable mirrors pumped by laser light with fre- 
quency u>o ■ Action of gravitational wave on such a system 
can be effectively described by means of effective forces 
acting on the mirrors and therefore changing dynamically 
the phase shift of outgoing light with respect to ingoing 
one. 

There exist two kinds of quantum fluctuations limiting 
the sensitivity of detector. They arc so called laser shot 
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FIG. 4: Schematic diagram of the quantum measurement de- 
vice. G is the classical observable {e.g. force) acting on the 
probe that is measured. Z is the output signal of the detector 
device, x is the linear observable {e.g. displacement) of the 
probe, and F is the linear observable of the detector which 
describes the back-action force on the probe. 



noise (SN) and radiation pressure noise (RPN). The first 
one originates from quantum fluctuations of electromag- 
netic wave phase which prevents from exact phase shift 
measurement and is, in essence, measurement accuracy. 
The second one, being a consequence of fluctuations of 
light amplitude, causes random radiation pressure force 
to move the mirrors and masks the measured signal force. 
As far as this noise is the direct consequence of measure- 
ment, it is also known as back action noise because of 
the back action of the measurement device (laser light) 
on measured quantity (phase shift). 

It is convenient to describe this system in terms of lin- 
ear quantum measurement theory formalism developed in 
[19] and thoroughly elaborated for use in gravitational- 
wave interferometers in [6]. Following this formalism we 
can represent our meter as two linearly coupled systems, 
probe (test masses) and detector (laser light and pho- 
todetectors). The schematic drawing of this equivalent 
linear system is presented in Fig. 4. Here x stands for 
some measured observable of the probe (mirrors relative 
displacement in our case), F is some observable of the 
detector through which it is coupled to the probe (ra- 
diation pressure force in our case), G is classical signal 
force being detected (GW action on the detector) and Z 
is the measured observable of the detector (output light 
quadrature in our case). Following [12] we write down 
the Hamiltonian of our system as: 

H = [{H v - xG) + H v ] -xF = F (0) + V , (4) 

where H(°> = [H-p — xG) + Ht> is considered as zcroth 
order Hamiltonian for both detector (marked by T>) and 
probe (marked by V) subsystems, and linear coupling 
between them V = —xF considered as a perturbative 
Hamiltonian. Using this approach, one can write down 
all of the observables of the system as a sum of unper- 
turbed zeroth order terms (marked by subscript (°') and 
perturbations (marked by subscript W) (see Eqs. (2.12- 
2.14) of [12]). In frequency domain these observables are 
read as: 



R FF (n)x^{n), 



(5a) 
(5b) 



cW(ft) = iW(n) + Lh{ti)/2 + i?^(s2)F (1) (s2) . (5c) 



Here Lh{fl) = R xx G{fl) is the GW signal proportional 
to metrics variation spectrum h(Q). Quantities Rab(Q) 
are frequency-dependent susceptibilities, in particular, 
Rzf (^) is the optomechanical coupling factor, 



R XX {Q) 



1 



mq 2 

is mechanical susceptibility of the SRI and 
RMO)- ra _ n a_2i 7 fi' 



(6) 



(7) 



is the optical rigidity. Z^ corresponds to shot noise of 
the laser light, x^ = Xt cc h stands for any displacement 
noise sources associated with the test mass reflecting sur- 
face with respect to its center of mass position, namely 
thcrmoelastic and Brownian noise of the mirror coatings 
and substrate, and 



F(°)(0) = F RPN (Q)+F tech (f2) : 



(8) 



where -Frpn(^) is the radiation pressure noise and 
-Ftech(^) is describing all the classical force noises, most 
notably suspension thermal, gravity gradient and seismic 
noises. 

Using these definitions one can now write down the 
output of the SRI reduced to metrics variation h units 
as: 

fcout(fi) = h{fl) + !{iMn)[FBPN(0) + ^tcch(ft)] 

+ [R xx {fl)R FF + lp SN (n) +X t ech(0)]} , (9) 

where 



xsn(^) 



R ZF {n) 



(10) 



is the normalized shot noise. 

Accounting for these definitions, one can write down 
now spectral density of the interferometer output noise 
as: 

s h {n) = s'^ nt {n) + s'; cch {n). (11) 

Here 

Ct( fi ) = ^(\R xx (m FF {n) + i|X SN (fi) 

+25R{[i? ra (r!)i? FF (fi) + l]*^ F (fi)}+|i? ra (fi)| 2 ^ PN (f7) 



(12) 



is the sum quantum noise spectral density, 

h n i + 2r 2 n 2 cos2/3 + r i 



2MJ 1 i 1 r 2 cos 2 (p + </>) + n 2 cos 2 ( 
is the shot noise (10) spectral density, 

2hMj 1 {r 2 + n 2 ) 



4 + 2r 2 ft 2 cos2/3 + r 4 ' 



(13) 



(14) 
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is radiation pressure noise .Frpn(^) spectral density, 



equal to: 



o fn\ . T sm{(3 + (f>) sin (f> 

SxF W = n rcos(p + 4>) + mco a <t> (15) 

is their spectral cross-correlation function, 

sum = ^(\R xx mR FF m + ii 2 sr h (^) 

+ \R xx (n)\ 2 s t F cch (n)) (16) 

is the sum technical noise spectral density, and ch (0) 
and Sl ech (fl) are spectral densities of non-quantum noise 
sources F tec h(Q) and cctech(^)- Optical losses influence, 
as shown in [14], can be accounted for by introducing 
effective quantum efficiency 77 of the readout photodetec- 
tor, which appears in Eq. (13). 

III. BINARY SOURCES 

The most popular and easy to implement criteria used 
to determine the optimal regime of GW detectors re- 
lates to evaluation of detection range for inspiraling bi- 
nary systems of compact objects such as neutron stars 
and/or black holes. This method is based on estimation 
of Signal-to-Noise-Ratio (SNR) using well known analyt- 
ical expression for spectral density of GWs emitted by 
system of two gravitationally bounded inspiraling astro- 
physical objects (See Sec. 3.1.3 of [20]): 



IM/)I 2 = 



G 5/3 n M 5/3 @ (/max _ f) 



12 



7/3 



(17) 



where M. = ^ 3 / 5 Af 2 / 5 is the so called "chirp mass" 
of the binary system constructed from reduced mass 
(j, = M1M2/M and total mass M = Mi + M 2 of the bi- 
nary system with components masses Mi and M2 corre- 
spondingly. One can readily see the indicative frequency 
dependence \h(f)\ 2 oc /~ 7 / 3 and inverse dependence on 
distance to the system squared r 2 . The upper cutoff fre- 
quency /max corresponds to the period of binary system 
rotation on innermost stable circular orbit (ISCO) when 
the system goes from quasistationary rotation phase to 
non-stationary merger phase. This frequency can be es- 
timated as: 



fn 



4400 x (Mq/M) Hz . 



(18) 



Given GW signal shape (17) and the noise spectral 
density (11), it is possible to write down the optimal 
SNR p which can be obtained on a given detector. As 
demonstrated by E. Flanagan and S. Hughes [21], SNR 
averaged over all mutual orientations between the detec- 
tor and source and over both polarizations of GWs is 



/™» | fe(/) |2 

J fm . m S h (2nf) 



df 



df 



2 g 5 / 3 A4 5 / 3 r f »- 

15tt 4 / 3 c 3 r 2 J fmja JV^SH^J) ' 



(19) 



where f m i n is the lower cutoff frequency at which binary 
system motion cannot be considered as stationary. In 
our calculations we will take / m i n — 10 Hz. 

In order to estimate detection range r one should set 
a threshold SNR po which defines the level of confidence 
in detection of GWs from binary system. Then detection 
range can be written as 



2 G 5 / 3 A4 5 / 3 



15 7T 4 / 3 C 3 



Pi 



/max 



df 



1/2 



fT/3 S \f) 



(20) 



Of course, all formulae we use here are obtained in 
the lowest Post-Newtonian order of general relativity 
[22] that definitely limits their application area to stellar 
masses systems and asymptotically flat space-time back- 
ground. However, for our purposes it is enough and the 
most significant for us feature of the above expressions is 
their relative simplicity. 

We calculated the detection range for standard (M = 
2.8-Mq) neutron stars binary system numerically con- 
sidering it as a function of 3 parameters: r G 
[500, 125001s" 1 , 13 G [-7T/2, tt/2], and (j> G [-tt/2, tt/2]. 
The number of values for each of the parameters was 193, 
giving in total 193 3 w 7.2 • 10 6 optical configurations. 
Distribution of points in parameter space was taken uni- 
form over angle variables (3 and </>, and logarithmic over 
variable F. For other parameters, we used the values 
planned for Advanced LIGO, see Table I. To account for 
technical noises, we used the noise budget also planned 
for Advanced LIGO and generated by BENCH software 
[23], see in Fig. 5. 

For each pair of r, (3 we have maximized SNR with re- 
spect to <fi, thus obtaining a function of only two param- 
eters -T, (3. The result of this calculation is presented in 
Fig. 6 as contour plots of normalized SNR /9nsns / Pnsns ■ 
Contours act in this plot as margins for regions in pa- 
rameter space where SNR is higher than the certain per- 
centage of maximal SNR Pnsns> b em g indicated in plot 
by point marked as "MAX". The parameters values for 
this maximal sensitivity point are listed in Table II. 

It is easy to note flat behaviour of SNR within a spa- 
cious range of parameters _T, (3. It arises due to two rea- 
sons. The first one is technical noises. It can be shown 
that in the absence of them, it is possible, in principle, 
to obtain arbitrary high values of SNR using deep and 
narrow well in quantum noise spectral density created 
by means of the second-order-pole regime of the opti- 
cal rigidity [9, 14], which corresponds to S = (J/4) 1 / 3 
and 7^0. Technical noise which has flat and smooth 
spectral dependence makes such excesses in the quantum 
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FIG. 5: Main classical noises planned for Advanced LIGO 
interferometer. 
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FIG. 6: Contour plot of sensitivity for standard NS-NS bi- 
naries Pnsns/Pnsns as a function of F, /?. Point "MAX" cor- 
responds to the sensitivity maximum, points "A" — 'D" - 
to typical suboptimal tunings shown in Fig. 7, point "O" 
- to double (NSNS+bursts) optimal regime, and points 
"P" , "Q" , "R" — to triple (NSNS+bursts+pulsar) suboptimal 
regimes. 



noise useless. Moreover, they increase the quantum noise 
at other frequencies and thus decrease the sensitivity. 

On the other hand, the integral character of criterion 
(20) allows significant variations of values of r and j3. In- 
crease of r decreases quantum noise level at low frequen- 
cies (/ < 100 Hz), but increases it at medium frequencies 
(100 Hz < / < 1000 Hz), and vice versa. Moderate pos- 
itive values of detuning S (and hence positive values of 
P) create resonance "wells" in the noise spectral density, 
but increase the low-frequency noise. As a result, the 
integral (20) does not change significantly. 

This consideration can be illustrated by Fig. 7, where 



FIG. 7: Quantum noise spectral densities optimized for stan- 
dard NSNS binary sources (point "MAX" in Fig. 6) and four 
typical sub-optimal spectral densities (points "A" — "D" (see 
parameters in Table. Ill) in Fig. 6). 





NSNS 


Bursts 


NSNS+Bursts+Pulsar 




MAX 


MAX 


"O" 


"P" 


"Q" 


"R" 


r 


2700 


2900 


3100 


4600 


6400 


6600 


P 


1.10 


0.57 


0.80 


0.52 


1.19 


1.39 


<t> 


-1.00 


-0.23 


-0.44 


-0.23 


-0.85 


-1.10 


p 2 


0.84 


0.74 


0.79 


0.80 


0.94 


0.97 


4>SRC 


1.48 


1.52 


1.51 


1.54 


1.53 


1.53 


7 = r cos j3 


1200 


2400 


2200 


4000 


2400 


1200 


8 = r sin 13 


2400 


1600 


2200 


2300 


5900 


6500 


Pns/pns x 


1.0 


0.988 


0.995 


0.989 


0.91 


0.84 


Pburst/ Pburst 


0.979 


1.0 


0.995 


0.989 


0.84 


0.75 


/ max 
Ppuls/ Ppuls 


0.21 


0.26 


0.25 


0.30 


0.40 


0.50 



TABLE II: Optimal ("MAX"), double-optimal ("O"), and 
triple-suboptimal ( "P" , "Q" , "R" ) parameters values for stan- 
dard NSNS binary, burst and typical periodic (pulsar J0034- 
0534) sources. 



quantum noise spectral densities for the optimal param- 
eters set (point "MAX" in Fig. 6) as well as for 4 typical 
sub-optimal ones (points "A" — "D") are plotted. 



IV. GRAVITATIONAL- WAVES BURSTS 

The next type of possible GW sources are supernovac 
explosions and stellar cores collapses, compact binary 
systems mergers [24] and other sources with not well 
modeled properties which are usually called GW bursts 
[15, 16, 17]. For these sources, we use the simple model 
of logarithmic-flat signal spectrum over the range of fre- 
quencies from /; to fh- By logarithmic- flat spectrum we 
mean that spectrum of GW signal h(f) is proportional 
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NSNS 


Bursts 




"A" "B" "C" "D" 


"A" "B" "C" "D" 


r 

a 
<$> 
p 2 

<j>SRC 


3161 7559 1487 1700 
-0.47 0.75 1.41 0.33 
0.97 -0.49 1.49 0.18 
0.85 0.97 0.73 0.75 
1.55 1.54 0.09 0.06 


3057 6182 2010 1438 
-0.70 0.49 1.24 0.20 
0.47 -0.29 -0.57 0.03 
0.84 0.95 0.77 0.73 
1.55 1.54 0.03 0.10 



TABLE III: Optical parameters for the SRI having quantum 
noise plotted in Fig. 6 and Fig. 9. 




FIG. 8: Contour plot of bursts sensitivity ptmrst /pSifrst as 
a function of T, (3. Point "MAX" corresponds to the sen- 
sitivity maximum, points "A" — "D" (see parameters in Ta- 
ble. Ill) — to typical suboptimal tunings shown in Fig. 9, 
point "O" — to double (NSNS+bursts) optimal regime, and 
points "P","Q","R" — to triple (NSNS+bursts+pulsar) sub- 
optimal regimes. 



to / x / 2 , that corresponds to constant numerator in ex- 
pression for SNR if integration is performed with respect 
to log/: 



2 

Pburst 



, fh \Kf)\ 2 df 
Z ' fl SWT 



log Sh dlQg f 

log/! sh^T) 



This way of defining SNR for burst events seems reason- 
able to characterize astrophysical signals with unknown 
spectrum structure so that contributions from frequency 
(time) ranges of different order are equal (for example, 
contributions from 10^ 100 Hz and 100^-1000 Hz should 
be the same). 

In Fig. 8, function pburst {r, (3) , calculated using the 
same algorithm as in the previous (NSNS) case, and nor- 
malized by its maximal value , is shown. Similar to 
the previous case, in Fig. 9 quantum noise spectral den- 
sities for the optimal parameters set (point "MAX" in 
Fig. 8 and for 4 sub-optimal ones (points "A" — "D") are 
plotted, and parameters values for the point "MAX" are 
listed in Table II. 



Total technical 
Burst optimal 
Burst suboptimals 




FIG. 9: Quantum noise spectral densities optimized for burst 
sources (point "MAX" in Fig. 8) and four typical sub-optimal 
spectral densities (points "A" — "D" in Fig. 8 with parameters 
given in Table. III). 



It follows from Fig. 8, that function pburst {r, (3) also 
has flat behavior within almost the same range of r, (3 
as Pns (r, (3) ■ The main difference between the NSNS 
and burst cases, invisible in Figs. 6 and 8, stems from the 
existence of cut-off frequency (18) and from more steep 
frequency dependence of the NSNS signal. As a result, 
the NSNS optimization procedure leads to smaller val- 
ues of angle (f> + f3, that reduces quantum noise at low 
frequencies, while the optimization with respect to GW 
bursts requires smaller values of <p that reduces quantum 
noise at high frequencies (compare Figs. 7, 9 and the cor- 
responding columns in Table 7). 

However, this difference is quite subtle, and it is evi- 
dent that regimes have to exist which provide good sen- 
sitivity for both these types of GW sources simultane- 
ously. In order to find them, we calculate values of r, (3, 
(j), which maximize the combined normalized sensitivity 



GNS+burst(A) — A 



PNS 
.max 
PNS 



2 + (1 - A) f £*E^V 



VPb, 



- . 

rst / 



(21) 



where A is a Lagrange multiplier which varies in the range 
[0,1]. 

The result is shown in Fig. 10, where parametric plot 
of Pbur S t(A)/Pb^t against PnsW/Pns* is presented. The 
leftmost point on this plot corresponds to A = 1, and the 
rightmost one — to A = 0. It follows from this calcu- 
lation, that the tuning exists where values of Pns and 
Pburst decrease both only by « 0.5% compared to their 
maximal values: 



Pns 



Pburst 



n max n max 
PNS Pburst 



0.995. 



(22) 



The corresponding values of parameters -T, (3, and <j> are 
listed in the column "O" of Table II. This point on the 
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FIG. 10: Parametric plot of Pburst(A)/p™,f r * t against 
Pns(A)/pns x 



HIGH-FREQUENCY PERIODICAL 
SOURCES OF GWS 



High frequency periodical sources of GWs, namely mil- 
lisecond pulsars, can be treated as very narrow-band al- 
most monochromatic sources with well defined central 
frequency 2/ pu i s [25]. For these sources, the detection 
range and the SNR are simply proportional to inverted 
square root of the noise spectral density at given fre- 
quency 2/ puls , 



'"puis Ppuls 



v/S /l (47r/ puls ) 



(23) 



Direct optimization of quantum noise in this case gives 
spectral density with very narrow and deep minimum at 
frequency 2/ pu i s , which is evidently non-optimal for the 
NSNS and burst sources considered above. Moreover, 
technical noise makes it useless to have very deep minima 
in quantum noise spectral density, limiting p pu i s by the 
value of 



max 
Ppuls (X 




FIG. 11: Quantum noise spectral densities optimized for 
NSNS sources (point "MAX" in Fig. 6), bursts sources (point 
"MAX" in Fig. 8) and for both of them (point "O" in 
Figs. 6,8). 



plane {-T, /?}, marked as "O", is shown also in plots 6, 8. 

In Fig. 11 the plot of quantum noise spectral density for 
this regime is presented together with plots of previously 
calculated optimal spectral densities for NSNS and burst 
sources. Note that all three curves are virtually indistin- 
guishable at low and medium frequencies (/ < 600 Hz), 
and the ones for bursts and combined NSNS+bursts are 
almost the same over all frequency ranges of interest. 



^cch(^/puls) 



(24) 



Therefore, we optimize the following "triple-purpose" 
function: 



GoNS+b urst-f puis (A,m) = 



i 



GNS+b urst (A) 



^puls 
/?puls 



(25) 

where < A < 1 and fi > are Lagrange multipliers. 
We took pulsar J0034-0534 [26] as an example of mil- 
lisecond pulsars, presumably emitting narrow-band high 
frequency GWs. Its barycentric rotational frequency is 
equal to fo rs 532.7 Hz and frequency of emitted GWs 
should be then /qw = 2/o ~ 1065.4 Hz. Among the 
high-frequency pulsars, the distance to this one is sig- 
nificantly smaller compared to other ones (0.98 kpc) and 
therefore is one of the most probable candidates for GW 
detection. 

The calculations results are presented in Fig. 10 as 
parametric plots of pns^Vpns* and Pbm-st(^) / Pburat 
against Ppu\s(^)/ p^is- It follows from these plots, that 
despite of two-dimensional character (two Lagrange pa- 
rameters A and n) of the optimization procedure, the 
results are virtually one-dimensional, because only very 
small trade-off between values of Pns an d Pburst is possi- 
ble (lines in both planes of Fig. 12 almost coincide). 

In Fig. 13 quantum noises for four characteristic 
regimes with optimal GNS+bm-st+puis are plotted. These 
plots demonstate how the optimization algorithm in- 
creases the sensitivity at the given frequency 2/ pu i s . The 
starting point "0" corresponds to obtained in Sec. V pa- 
rameters set, optimized for NSNS and bursts signals, 
see Figs. 6, 8. First, the optimization algorithm tries to 
suppress the quantum noise in all high-frequency range 
(/ > 10 3 Hz) by increasing r and decreasing /3, see point 
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FIG. 12: Parametric plots of Pns(/x)/Pns x (top) and 
Pburst(/i)/pwst (bottom) against p pu is(A)/p™ u a 1 g for two opti- 
mizations regimes: with priority to pns (A = 1) and to pburst 
(A = 0). 



"P" in these plots. At this stage, noticeable gain in pul- 
sars sensitivity (~ 1.5) can be obtained with negligibly 
small (~ 1%) sensitivity loss for NSNS and bursts sources 
(see Table II). Then, the optimization algorithm starts 
to "grow" local minimum at frequency 2/ pu i s , by increas- 
ing back P in such a way that S = _Tsin (3 — > 47r/ pu i s , see 
points "Q" and "R" . At this stage pns and pburst start to 
decrease noticeably (by tens of percents). However, p pu i s 
increases several times at this stage. 



VI. CONCLUSION 

The results of this paper rely heavily on the esti- 
mates of the technical noise predicted for Advanced 
LIGO. These estimates almost definitely will be subject 
to changes during the next few years, however, it is im- 
probable that technical noise estimates will change sig- 



nificantly. Therefore, all specific values obtained here 
should not be considered as final ones. 

The main result of this paper is not these values, 
but the conclusion that regimes of the signal-recycled 
interferometer exist which can provide good sensitiv- 
ity for both binary inspiraling and burst gravitational 
wave sources. Moreover, "triple-purposes" regimes are 
also possible, which provide significant sensitivity gain 
for high-frequency periodical sources (millisecond pulsar) 
with only minor sensitivity degradation for binary inspi- 
raling and bursts GW sources. 

Calculations presented in this paper show that in order 
to obtain good sensitivity for binary, burst and, to some 
extent, high-frequency periodic sources, it is necessary to 
use large values of interferometer bandwidth 7 ~ (2-^4) x 
10 3 s _1 3> 2 x 71TOO s _1 with significant positive detunings 
if~2x 10 3 s" 1 , see Table II. These tunings give smooth 
broadband shape of quantum noise curves, dictated by 
the technical noises, especially by the mirrors thermal 
noise which has very flat frequency dependence in GW 
signals spectral range. 
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